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Abstract 

In this paper, we study a nonlinear two-dimensional two-particle 
quantum tunnel effect with dissipation in diatomic H-H system. We 
use an instanton technique based on the path integral approach, and 
present analytical solutions, in the approximation of linearized reac- 
tion coordinates and harmonic potential. It appears that the tunnel 
effect leads to a specific correlation between two electrons interacting 
with each other by Coulomb potential. This effect can be viewed as a 
mechanism supporting the isoelectronium state recently proposed by 
Santilli and Shillady within the framework of isochemical model of the 
hydrogen molecule. We show that, under a condition of effective two- 
dimensional motion, contribution of the quantum tunnel correlated 
configurations is essential in the electron dynamics. Temperature de- 
pendence of the probability of the two-electron tunnel transitions per 
unit time has been studied. 
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1 Introduction 



In a recent paper, Santilli and Shillady U suggested a new isochemical model 
of the hydrogen molecule characterized by a bond at short distances of the 
two valence electrons into a singlet quasi-particle state called isoelectronium. 
Namely, the isochemical model suggests introducing of attractive short-range 
Hulten potential interaction between the electrons, in addition to the usual 
Coulomb repulsive potential. An origin of the short-range attractive potential 
has been assumed due to a deep overlapping of wavepackets of the electrons 
giving rise to nonlinear and other effects. 

Both the cases of unstable and stable isoelectronium corresponding to 
the general four-body and the restricted three-body problems, respectively, 
have been analyzed, in Born-Oppenheimer approximation. Standard Boys- 
Reeves numerical three-body and four-body calculations Jl|, exact analytic 
three-body solution f2|, and analytic four-body Ritz variational study || 
of the model have been carried out. The Hulten potential Vh, which con- 
tains two parameters, appeared to lead to rather complicated calculations 
so that the simplified potentials, the Gaussian-screened-Coulomb potential 
V g and the exponential-screened-Coulomb potential V e , have been used in 
the above papers, to mimic the Hulten potential at long and short distance 
asymptotics. A comparison of the obtained results with experimental data 
on H 2 molecule showed that the approximation of stable pointlike isoelectro- 
nium (the three-body problem) does not fit the data while for the unstable 
isoelectronium case (the general four-body problem) exact representation of 
the experimental binding energy and bond length of H 2 molecule can be 
achieved by fitting the correlation length parameter of the potential. The 
ground state energy based estimation || showed that the contribution of the 
attractive interelectron potential V e (i.e., the contribution of the potential 
related to the isoelectronium) should be about 1%...6% of that of the usual 
repulsive interelectron Coulomb potential, to meet the experimental data. 
In ref. pi certain relation between the two parameters of the Hulten poten- 
tial was established, thus only one parameter has been used for the fitting. 
This relation arose from consideration of the two-electron system governed 
by the Hulten potential, with the assumption that the electron pair could 
form one-level bound state with zero ground state energy. 

The effect of electron correlations made by the introduced attractive 
short-range interelectron potentials (i.e., the unstable isoelectronium), with 
the correlation length parameter value of about 0.01 bohr appeared to be 



1 



helpful in achieving exact representation of the experimental data on H 2 
molecule. Also, within the Boys-Reeves framework the isoelectronium based 
approach appeared to provide an essentially (at least 1000 times) faster com- 
puter calculations |IJ, in comparison to a standard C.I. calculation. This 
approach has been extended to the case of other diatomic molecules |jj and 
water molecule 

In the present paper, we study two-electron correlations arising, alter- 
natively, from the two-dimensional tunnel effect with dissipation which can 
take place in diatomic molecules subjected to strong external electromag- 
netic fields. More specifically, we take H2 molecule as a simple example and 
assume that the three-dimensional potential of the model can be effectively 
reduced to a two-dimensional one under the action of external fields. Also, 
we assume that the H2 molecule is a part of molecular association. Due to 
obtained results, the two-particle tunnel effect leads to a contribution which 
resembles that of the above isoelectronium correlations in two dimensions. 
However, we should to emphasize that we do not introduce any attractive 
interelectron potential. Namely, the characteristic two-electron correlation 
arises naturally as a consequence of the two-dimensional two-particle tunnel 
effect with dissipation. 

Present consideration is strongly motivated by recent studies of a light gas 
called magnegas™, consisting due to a standard chemical analysis mainly 
of carbon monoxide CO (41%) and hydrogen H 2 (48%), produced by San- 
tilli's PlasmaArcFlow™ chemical reactor |3, H|. Gas-chromatography mass- 
spectrometry and infrared spectroscopy data obtained at room temperatures 
clearly indicate presence of high molecular mass species (up to 1000 a.m.u.) 
in magnegas, in a macroscopic percentage. They have no strong infrared sig- 
nature of conventional chemical bonds, except for that of carbon monoxide 
CO and carbon dioxide CO2, which is present in magnegas in small percent- 
age. The GC-MS/IR search results using a library of about 138,000 chemical 
species did not indicate any matches with these high molecular mass species. 
The main hypothesis on the origin of these species called magnecules is that 
they are formed from the usual molecules (CO, H 2 , etc.) bonded to each 
other in some way, under the initial action of strong external electromag- 
netic fields produced by underwater arc in the PlasmaArcFlow reactor pf. 
We refer the interested reader to refs. H |(J for more detailed information on 
magnegas and PlasmaArcFlow reactor, and their applications. 

In the present paper, we conjecture that the two-dimensional tunnel ef- 
fects play an essential role at the stage of forming of the above mentioned high 
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molecular mass species (magnecules) inside the reactor, where the influence 
of strong external electromagnetic fields takes place. Physical and chemical 
processes inside the reactor concern specific plasmachemical studies and can 
be elaborated elsewhere. Here, we note only that the importance of studying 
of the tunnel effects is grounded on the fact that a tunnel effect is known to 
be as a mechanism in forming of formaldehyde (CH 2 =0) polymer structure. 

We concentrate on the study of consequences implied by the eventual 
'magnetic freezing' (lowering of the effective temperature due to the effect of 
external magnetic fields) of the hydrogen molecule in molecular association. 
The magnetic freezing is characterized by an effective reduction of the three- 
dimensional potential to some two-dimensional one ('polarizing of orbits'). 
We should to emphasize that the system under study is not a usual H 2 
molecule since we assume the effective two-dimensional treatment of electron 
dynamics and the presence of neighbor molecules (molecular association) 
serving as a heat bath; herebelow, the system under study is referred to as 
H-H system. 

We use oscillator potential as an approximation to superposition of Cou- 
lomb potentials in the H-H system in the molecular association. Such an 
approximation is naturally used in low-temperature chemical kinetics, with 
parameters of the oscillator potential being adjusted to reproduce a ground 
state energy level. Also, we account for the interelectron interaction by 
using the Coulomb potential. We note that the ground state energy of the 
"magnetically frozen" H 2 molecule (i.e., the H-H system) is lower than that 
of the ordinary H 2 molecule. As the result the binding energy of the H-H 
system is bigger than that of H 2 molecule. Detailed study of the effect of 
external magnetic fields is out of the scope of the present paper and will be 
made elsewhere. 

Set up of the model is the following. We consider two-dimensional (pla- 
nar) motion of two electrons around fixed nuclei of two H atoms. The elec- 
trons thus move in a nonlinear potential U(q\,q2), where the reaction co- 
ordinates qi and g 2 of the electrons in the H-H system can be linearized in 
the first order approximation as shown in Fig. 1. The linearization assumes 
introducing of four centers, instead of the original two centers (protons). 

The resulting approximate nonlinear potential U(qi,q2) represents two- 
dimensional potential surface modelling the two-dimensional potential of the 
original system (two protons and two electrons moving in the plane). This 
surface has the highest points at a one-half of the proton-proton distance, 
and four lowest points at the proton positions. 
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Figure 1: The linearized reaction coordinates qi and q2- 



Dynamics of two interacting electrons in this two-dimensional potential is 
due to the ordinary quantum mechanics to which the energy barrier between 
the protons makes a specific contribution: two-dimensional quantum tunnel 
effect of two-electron transitions. The tunnel component of the probability 
of two-electron transitions becomes essential in the case when the ordinary 
overbarrier two-electron transitions are suppressed. 

Clearly, the tunnel transitions can greatly affect both the character and 
strength of the chemical bond in diatomic molecule. As a consequence, they 
can give an essential contribution to electron dynamics in the molecular as- 
sociation consisting of such molecules. 

We remark that usually a tunnel effect is considered for the case of one 
particle in one dimension, which is a simplest case studied to much extent. 
In this paper, we consider the tunnel effect with dissipation for the case of 
two interacting particles in a two-dimensional potential. 

The tunnel correlation of electrons is considered on the background of a 
formed two-dimensional potential of the H-H system which is assumed to be 
a part of the molecular association. Therefore we study temperature effects 
by introducing interaction of the H-H system with a heat bath. 



2 The two-dimensional tunnel model 

Applications of the theory of quantum tunnelling with dissipation [|7|- ||14| in 
chemical kinetics can be found in the literature [15fl-|17l. Below, we turn to 
determining of the two-dimensional nonlinear potential for two interacting 
electrons. 

For the case of non-interacting electrons the potential energies of first 
and second electrons, as functions of the reaction coordinates, q\ and q2, are 
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taken in the following form: 

U(qi) = \uj\ qi + a) 2 6{- qi ) + l -uj\ qi - a) 2 9( qi ), 

U(q 2 ) = \^\q 2 - af9{q 2 ) + l -u\ qi + a) 2 9(-q 2 ), (1) 

where 2a is the distance between two nuclei in the H-H system (e.g., 2a ~ 1.40 
bohr, for H 2 molecule in vacuum), mass of electron is m = 1, and u is a 
frequency. In Fig. 2, the potentials U(qi) and U(q 2 ) are shown. 
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Figure 2: The potentials U(q\) and U(q 2 ) as functions of linearized coordi- 
nates. 



The interaction between two electrons is taken in the following form: 

V int (qi, q 2 ) = ~&{<li - q 2 ) 2 , (2) 

where a is a positive parameter, and (gi — <j 2 ) is distance between the elec- 
trons. Note that Vi nt has the form of attractive harmonic potential. Indeed, 
in the xy-plane the x-component of the mutual distance of two interacting 
electrons in two dimensions are much bigger than the y-component, % ^> a, 
as shown in Fig. 3 

Thus, the interaction potential can be expanded in small parameter (g ly — 
Q2y) 2 /q 2 , where q iy and q 2y are coordinates of tunneling. The Coulomb re- 
pulsion potential between the electrons, with electric charges — e, is then 



V c = 
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Figure 3: Transition of two interacting electrons. 



where Eq is the dielectric constant. The negative (second) term in the series 
expansion (|3[) is the effective attractive potential (Q) so that the constant a 
is found as 

e 2 

a = (4) 

The first term, e 2 /(eo ( ?o), is responsible for repulsion between the electrons, 
and is constant along the reaction coordinates so that it can be absorbed by 
redefining U(qi) and Ufa). Obviously, the interaction potential Vc is always 
repulsive but one can see that only its effective attractive part given by Eq. 
(fj), as a function of the reaction coordinates, appears to contribute to the 
reaction dynamics. 

The following remark is in order. An attractive character of the inter- 
electron potential V int (qi,q2) resembles the attractive interelectron Hulten 
potential V^(r) of the isochemical model 0. The difference is that Vi nt 
is formulated in terms of reaction coordinates, has a harmonic form, and 
naturally arises as an effective long-range potential in the two-dimensional 
approach. So, Vi nt can be viewed as a potential essentially supporting the 
isoelectronium state in two dimensions. 

Thus, the general form of the total potential is a sum of U(qi), U(qi), 
and Vi„ 4 (gi,g 2 ), 

U (qi, q 2 ) = U ( gi ) + U (q 2 ) + V int ( qi , q 2 ), 
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which we rewrite as 



U{qi,q 2 ) 



( qi+ aye(-q 1 )e(q 1 ) + (q 2 -aye(q 2 )e(-q 2 )- Y ( qi -q 2 ] 

(5) 



Here, we have denoted 
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which is a dimensionless parameter, a* < 1. In Fig. 4, a section of the 
potential (|5|), at some energy level, is schematically depicted. In Figs. 5 
and 6, three-dimensional plots of the potential, at two different values of 
the renormalized interelectron coupling parameter, a* = 0.1 (weak coupling) 
and a* = 0.5 (strong coupling), are shown. Here, the ovals corresponding 
to sections of the potential by a horizontal plane at some energy level are 
depicted for the reader convenience. We recall that q\ and q 2 are the reaction 
coordinates (not coordinates of real configuration space of the system) . 
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Figure 4: Level lines (1,2,3,4) of the potential U(qi, q 2 ); a* 
coupling parameter. 



is the renormalized 
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Figure 5: A three-dimensional plot of the potential U(qi,q 2 ); a* = 0.1, a = 
0.7. 



3 Probability of the two-electron transition 

Since the hydrogen molecule is assumed to be a part of an association of 
molecules (magnecule), we introduce interaction of the H-H system with 
external oscillators, a heat bath. Dynamics of the heat bath is defined by 
the oscillator Hamiltonian 

H hb = lUP* + ^Q*)- (6) 

i 

We assume that each of the electrons linearly interacts with the oscillators, 
namely, the interaction potentials are 

V e %(qi, Qi) = qi £ QQi, V^ hb (q h Q t ) = q 2 £ QQ,, (7) 

i i 

The probability of the two-electron transition per unit time is given by 

ImZ . , 

= ReZ' < 8) 

where, as usual, the most important part is the exponential part of T. At 
zero temperature, for the case of metastable levels we have 

r = -2Im£, E = E -^T, (9) 



Figure 6: A three-dimensional plot of the potential U(qi,q 2 ); a* = 0.5, a = 
0.7. 



from which Eq. (§|) follows as the generalization to the case of finite temper- 
ature T. Indeed, 

2Eexp(-E 0l /T)lmE i 2Tim£exp(-£ i /T) 

J2ex V (-E 0i /T) ReEexp(-E i /T) ReZ' K J 

i i 

Here, i runs over energy levels of the metastable states, and Z is the statistical 
sum of the total system. The imaginary part of Z corresponds to a decay of 
the energy levels. 

To calculate V it is convenient to represent Z as a path integral 0-[14] 
Z = Iij Dq 1 Dq 2 DQ i ex P [-S(q 1 ,q 2 ,Qi)], (11) 

i 

where S is the action of the total system. An imaginary part of Z corresponds 
to the decay of the initial energy levels. One can perform functional integral 
over Qi exactly, and obtain 



Z = Il/ D( li D( l2 exp[-S(q 1 ,q 2 )}, 



where 

13/2 



S(qi,q 2 )= J dr{^q 2 1 + ^q 2 2 + V(q 1 ,q 2 ))+ (12) 



0/2 

2 / ^(r-r')(gi(r) + g 2 (r))( gi (r') + g 2 (r'))}, 



0/2 

where 

1 oo 

W = 77 E ^K)^ (13) 

A 3 71= 



-oo 

is Green function for the oscillators, with 



^) = -E;i?3. (14) 



7" W i + 



/3 = h/ksT is the inverse temperature parameter, and is Matsubara's 
frequency. 

The quasiclassical trajectory minimizing the action S in two-dimensional 
space is defined by the equations of motion, 



0/2 

-qi+nlqi+a 1 q 2 + J dr K(T~T , )(q 1 (T , )+q 2 (T')) +u 2 ae(-q 1 )-u 2 a6(q 1 ) = 

-0/2 

(15) 



and 

0/2 



-g 2 +^g 2 +a 1 gi+ / rfr , K(r-r , )(g 1 (r , )+g 2 (r , ))-cJ 2 a^(g 2 )+w 2 a0(-g 2 ) = 0. 

(16) 



-0/2 



Here, the kernel K is defined by 

K{r) = - £ in^\ (17) 

P n=—oo 

where £ n is determined from the redefined Eq. flliD , namely, 

%) = -E^ + ^ (i g ) 

i ^ 

This redefinition has been made in order to extract a zero mode term. 

We seek for solutions of the set of Eqs. (|15|) and ([16]) in the form of 
Fourier series expansion in frequencies u n 
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1 00 1 
?i = 77 E tf^, * = /? E (19) 

P n=— oo P 



n=— oo 
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The renormalized frequency, fio an d renormalized coupling constant a% are 
defined by 

^o = ^ 2 -E^|-«' ai = a-Y,% (20) 



LOT ; LOt 



Inserting Eqs. fll9|) into the set of equations (|l~5l) and ([TBI) , we get the 
following equations: 

(i) , (2) 4^ 2 ae 



fig + ai 



Y (i) _ J2) _ 2a; 2 a/3 8u; 2 To 
fig — «i — a i 



= + (21) 

for the case n = 0, and 



(i) , (2) _ 4uTa(sin zyri - sin z/ n r 2 ; 



i/ n (i^ + fig + ax + 2£ 



ni 



(i) _ ( 2 ) _ 4u 2 a(sin v n T X + sin u n r 2 ) 

Qn Qn ~ ( 2 i n2 \ ' v / 

z/ n (z/ 2 + fig - ai) 
for the case n ^ 0. Here, we have denoted 

£ = Ti-t 2) r = ^(n + r 2 ), (23) 

and ri and r 2 are the values of time r at which first and second electron, 
respectively, passes through the maximum of the potential barrier. Namely, 
7~i and r 2 are determined from the following equations: 

? 1 (n) = 0, g 2 (r 2 ) = 0. (24) 

These two equations can be used to change arguments of the above functions 
9. As the result, the 9 functions become depending on the time values T\ 
and r 2 , instead of the coordinates q\ and g 2 , and the Eqs. (113 ) and flTED 
take a linear form. Note that, quasiclassically, the time values ±ti and ±r 2 
correspond to the times when first and second electron, respectively, are at 
the top of the potential barrier. 

Inserting the trajectories defined by Eqs. (|I9D , (FH), and fl22| ) into Eqs. 
(|T2|), we obtain the quasiclassical (instanton) action in the form 

8cu 4 a 3 r 4w 4 aV 16cu 4 a 2 r 2 



fig - oc x /3(fig + «0 /?(fig - ai) 
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32w 4 a 



4„2 oo 



E 

71=1 



sin 2 z/„r + cos 2 v n ej2 sin 2 ^ n e/2 + cos 2 u n r 



+ 



^ + fig - ai ) ^2 + fig + a x + 2£ n ) 



(25) 



The quasiclassical action S describes probability of the two-particle tunnel 
transition per unit time, with an exponential accuracy. For the case of an- 
tiparallel motion of the electrons, the above instanton action S, at £ n = 
(the heat bath effects are ignored here), takes the form 



„ , 1 1 i/-, In sinh(|£|\/I - a*) . , . . 
S = -2uja 2 < e 1 ) H - smh e 



1 - a* 



+ COSh(l ! l ^ ) + 1 ■ mShW ' ~ r) S^^ (26) 



(1 - a*) 3 / 2 
4 



sinh ( /?* \/l — a*) 
(coshe - l)(cosh(/3* - r) + cosh/3*) 1 



sinh p* 



Here, we have denoted 



a = -j, P = TT-, e = (ri - r 2 )u;, r = (n + r 2 )u;, 
and the parameters e and r obey the following set of equations: 
— sinhe[coth /5*+coshr coth (3*— sinh r ]+^ ~ sinh(e\/l — <**) [coth(/3* vT 



- cosh(r v / l - a*) coth(/3*Vl - a*) + sinh(r v / l - a*)] = 0, 
1 



- + (cosh e — 1) (sinh r coth (3* — cosh r) + cosh e (27) 



1 — a 

+ - 1 {[cosh(eVl - «*) + l][sinh(r v / l - a*) coth(p*Vl - a*) 

— cosh(e\/l — d*)] — cosh(£\/l — a*) \ . 
Solution of the Eq. (|7]) is 

£ = ( Tl - r 2 )u; = 0, a < ^u 2 , V/3, 



(28) 



7i = r 2 



2a; 4' 
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At low temperatures, u/3 3> 1, we then have 



a' 



e -rVT^^ a * {1 _ a y Ji_ ( i_ qT( __^__i)^ . ( 29) 



i (a* - e - r ^n 

g & r^j 

I - a* e -rVi-«* 

where 

_ yi - a* 
7 ~ 1 - VI - a*' 
The approximate solution (|29| ) is valid at 

1 2a _ 

and (3 > (3 C , where the critical temperature is 

1 , \ a*(l-a*) 7 

As one can see, there arises a characteristic critical temperature, /3 C = Tijik^T,}, 
of the system, and it depends mainly on the value of the interelectron cou- 
pling parameter a*. 

At e = 0, for the case of symmetric two-dimensional potential the solution 
(H) yields 

(1 - a*) 3 / 2 4 ^ ( } 

At £ 0, the resulting expression for the action 5 is cumbersome so that we 
do not represent it here. However, it can be shown that, at (3 > f3 c , 

S(e=0) < S( £ ^Q) (32) 

so that the value e = minimizes the action. As the result, the values e ^ 
correspond to splitted trajectories which appear to be unstable; see Fig. 7. 

The instanton solution for e = corresponds to a strongly correlated 
antiparallel motion of the two tunnelling electrons. Namely, the two electrons 
simultaneously pass the top points of the potentials. This correlation of 
the two electrons depends on the value of the electron-electron (attraction) 
coupling constant a* naturally arising within the framework of our model. 

Such a correlation resembles the isoelectronium correlation introduced in 
ref. ]I] where the correlation was introduced within a different approach and 
is governed by the short-range attractive electron-electron Hulten potential. 
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Figure 7: A schematic picture of the tunnel trajectories at e = (solid line) 
and at e ^ (dash lines). 

4 Conclusions 

At f3 > j3 c , the basic single trajectory is splitted into two trajectories as 
shown in Fig. 7 (dash lines). In contrast to the case of parallel transition, 
this splitting (bifurcation) takes place at any values of the parameters of 
the potential U(q 1: q 2 ). Also, at f3 > f3 C) we have S^ ) > 5 , ( £ =o) so that 
the trajectory with e = make a biggest contribution. At (5 < f3 c , the two 
trajectories become a single trajectory characterized by q\ = —q 2 . 

We see that in contrast to the case of one tunnelling particle at which only 
one trajectory (instanton) is realized, in the case of two tunnelling particles in 
the two-dimensional potential the situation is more complicated. As we have 
shown in general there are two types of trajectories, a single basic trajectory 
and splitted (degenerate) trajectory, both making contribution to the tunnel 
effect. The splitted trajectory is characterized by a nonsimultaneous {j\ ^ r 2 ) 
transition of the two particles through the top points of the barriers, so 
that the time-correlation between the particles' motion is lost. It is highly 
remarkable that in the case of symmetric two-dimensional potential such 
splitted trajectories in both the cases of parallel and antiparallel motion of 
the two particles are unstable. At small a* and (5 < f3 c , the degenerate 
trajectory is not realized so that there is only the basic single trajectory 
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(qi = —qi). This trajectory is characterized by a simultaneous (ti = r 2 ) 
transition of the two particles through the top points of the barriers, so that 
we observe a specific correlation between the electrons. 

The chosen form of the interelectron interaction does not affect the motion 
of the center of mass (gi = q 2 )- So, for the parallel transition along the basic 
trajectory (g x = q 2 ) the quasiclassical (instanton) action does not depend on 
the coupling parameter. In this case, 

S = Aua 2 tanh — . (33) 
4 

Since the motion with maximal value of the relative coordinate, i.e. q% = —q 2 , 
is energetically preferable it becomes clear why for the parallel transition 
along the degenerate trajectory the action decreases with the increase of the 
coupling parameter, while for the antiparallel transition along the degenerate 
trajectory the action increases with the increase of the coupling parameter. 

Thus, in the case of antiparallel tunnel transition of two interacting elec- 
trons the preferable trajectory is a single basic one which is characterized 
by strong (isoelectronium-like) correlation between the electrons (tj = r 2 ). 
Here, the uncorrelated motion ("decay" of the isoelectronium-like state) is 
suppressed because it makes bigger contribution to the instanton action. 

As to a heat bath, it appears to be possible to study analytically the effect 
of, for example, one local mode ul to the two-electron transition probability. 
The result is that the heat bath does not change the qualitative conclusions 
made above. Temperature dependence of the action S for the antiparallel 
motion case (e = 0) at two different values of the renormalized coupling 
parameter, a* = 0.1 and a* = 0.5, is shown in Fig. 8 

In this paper, we used the assumption that dynamics of the two inter- 
acting electrons can be taken approximately as a two-dimensional. This is 
due to the viewpoint Jl| that strong external magnetic field confines, to some 
degree, the usual three-dimensional motion of the electrons so that the con- 
sideration of two-dimensional effects in molecules, as well as in associations 
of molecules, becomes important. The main two-dimensional effect, namely, 
the two-electron tunnel effect with dissipation, for the antiparallel transition, 
has been studied in this paper. 

By using the approximations of one-instanton linearized reaction coordi- 
nates and an ideal gas of instanton-antiinstanton pairs ((3u) <C 1QU barrier / , 
we have shown that the correlated isoelectronium-like motion of the two tun- 
nelling electrons in two-dimensional H-H system with dissipation takes place, 
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Figure 8: The action (e = 0, antiparallel motion) as a function of inverse 
temperature (3 = h/k B T; (a) a* = 0.1, (b) a* = 0.5. 
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and appears to be the most stable configuration. Thus, the tunnel correla- 
tions can be viewed as a mechanism supporting the idea of isoelectronium 
states introduced by Santilli and Shillady 

Detailed study of the effect of strong external magnetic fields, as well as 
the implications of the results of this paper, will be made in a subsequent 
paper. 

Application of the presented instanton formalism to carbon monoxide 
molecule CO is of much interest as well. Also, study of the transition rates, 
T = e~ s , for the case of association of molecules (intramolecular tunnelling) 
could be important in investigating electron charge distribution and bonds 
in the association of molecules (magnecules) . 
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